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Physical notions of stochastic resonance for potential diffusions 
in periodically changing double-well potentials such as the spectral 
power amplification have proved to be defective. They are not ro¬ 
bust for the passage to their effective dynamics: continuous-time 
finite-state Markov chains describing the rough features of transi¬ 
tions between different domains of attraction of metastable points. In 
the framework of one-dimensional diffusions moving in periodically 
changing double-well potentials we design a new notion of stochastic 
resonance which refines Freidlin’s concept of quasi-periodic motion. It 
is based on exact exponential rates for the transition probabilities be¬ 
tween the domains of attraction which are robust with respect to the 
reduced Markov chains. The quality of periodic tuning is measured by 
the probability for transition during fixed time windows depending 
on a time scale parameter. Maximizing it in this parameter produces 
the stochastic resonance points. 

0. Introduction. One of the simplest and earliest stochastic climate mod¬ 
els goes back to [1] and [17]. It intends to give a qualitative explanation of 
glacial cycles and is based on a deterministic differential equation for the 
global mean temperature expressed through a balance between the albedo- 
driven absorbed and the black-body type emitted radiative energies. A pe¬ 
riodic exterior forcing comes from the slowly fluctuating solar constant and 
is due to Milankovich cycles caused by the gravitation of big planets. Only 
the addition of a stochastic term as a second forcing makes spontaneous 
transitions between the otherwise isolated metastable states of temperature 
possible. The resulting stochastically and periodically perturbed differen¬ 
tial equation was capable of describing at least one characteristic aspect of 
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experience: the typically short and abrupt transitions , observed before by 
Kramers [14] in reaction-diffusion phenomena. The model was soon strongly 
disputed. Despite its lack of realistic assumptions, the concepts underlying 
the model brought to light the phenomenon of stochastic resonance. Roughly 
speaking, a periodic (input) system subject to random perturbations is in 
stochastic resonance, if the noise intensity is tuned in such a way that the 
random periodic output is optimal. A very lively research field developed 
around this concept, drawing numerous examples from a wide spectrum of 
areas (see [9, 10] for a survey). 

The mathematically precise understanding of the phenomenon is still un¬ 
der discussion at the time this paper is written. The first approach is done 
in [7], where the deep large deviations’ theory of [8] is employed to produce 
a notion of stochastic resonance explaining the phenomenon in the small 
noise limit as approximating the periodic hopping between the energetically 
most favorable states in the landscape provided by a periodically weakly per¬ 
turbed potential with finitely many local minima. In this sense stochastic 
resonance can be understood as the ability of the system to undergo quasi- 
periodic motion in the limit of small noise intensity. Let us briefly recall 
this interpretation more precisely. If noise intensity is e, in the absence of 
periodic exterior forcing, the exponential order of times at which successive 
transitions between metasable states occur corresponds to the work to be 
done against the potential gradient to leave a well. This fact, heuristically 
derived by Kramers and Eyring ( Kramers ’ time), is shown with mathemat¬ 
ical accuracy in [8]. The attractor basins are subdivided into a hierarchy of 
cycles with main states corresponding to the deepest among the cycle states. 
In the presence of periodic forcing with period time scale e My/£ , in the limit 
e —> 0 transitions between (the main states of) cycles with critical hopping 
work close to p will be periodically observed. Transitions with smaller crit¬ 
ical work may happen, but are negligible in the limit. Those with larger 
critical work are forbidden. In the simplest case of two minima of potential 
depth V and v, v < V, the role of which switches periodically at time T, 
for T larger than e v ' £ the diffusion will be quasi-deterministic, that is, close 
to the deterministic periodic function following the location of the deepest 
well. 

Quasi-periodicity captures an important aspect of stochastic resonance, as 
it provides conditions under which stochastic trajectories are able to exhibit 
periodic behavior. Yet, physics literature (see [9, 10]) stipulates that stochas¬ 
tic resonance not only explains conditions for stochastically periodic behavior 
but comprises its optimality in a sense quite similar to the resonance notions 
of wave dynamics. In classical optics resonance is understood as the optimal 
amplitude of the response of the system to periodic excitation. In the same 
way, a stochastic resonance point is claimed to explain optimal periodic tun¬ 
ing of the stochastic trajectories of the diffusion responding to deterministic 
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periodic excitation. Amplitude as a measure of quality of periodic tuning is 
replaced by signal-to-noise ratio or spectral power amplification (see below). 
Numerical simulations as, for example, in [16] clearly support the optical 
evidence that beyond the threshold described by Freidlin [7] at which quasi- 
deterministic behavior becomes possible, for different noise intensities quite 
different qualities of periodicity of the random trajectories can be observed. 
There are parameter ranges for T in which random trajectories follow quite 
well the deterministic shapes of excitation curves. But as T gets even big¬ 
ger, many short excursions to the wrong well during one period may occur. 
They will not count on the exponential scale on which quasi-periodic mo¬ 
tion is measured, but trajectories will look less and less periodic. Physicists’ 
quality measures for tuning therefore cannot be explained on the basis of 
quasi-deterministic motion alone. 

The thesis by Pavlyukevich [13] and [18] presents an attempt to pro¬ 
vide a mathematically sound underpinning of physical notions of stochastic 
resonance based on optimality of periodic tuning—as opposed to the trajec- 
torial analysis of Berglund and Gentz [2] containing very fine estimates on 
relaxation times. The physical concepts are mostly based on comparisons of 
trajectories of the noisy system and the deterministic periodic curve describ¬ 
ing the location of the relevant metasable states, averaged with respect to 
the equilibrium measure. In the simple one-dimensional situation considered 
above the system switches between a double-well potential state U(x) with 
two wells of unequal depths V and v, v < V, during the first half period, and 
U(—x) for the second half period. The total period length is T, and stochas¬ 
tic perturbation comes from the coupling to a white noise of intensity s. 
The most important measures of quality studied are the spectral power am¬ 
plification, the related signal-to-noise ratio or the entropy of the equilibrium 
distribution. In particular, the first two mentioned play an eminent role in 
the physical literature. They mainly contain the L 2 average in equilibrium 
of the spectral component of the solution trajectories corresponding to the 
input period T, normalized in different ways. These measures of quality are 
functions of £ and T, and the problem of finding the resonance point, for 
example, consists in optimizing them in e for fixed (large) T. 

Let us briefly explain a striking shortcoming of these concepts of opti¬ 
mal periodic tuning which made us look for different ones. The first step to 
find optimal tuning intensities e(T) for large T consists in reducing the dy¬ 
namics of the diffusion to the interwell motion, that is, the pure transitions 
between the potential minima. In the physics literature, this corresponds to 
the reduction given, for example, by [15]. One ends up with continuous-time 
two-state Markov chains with transition probabilities corresponding to the 
inverses of the diffusions’ Kramers-Eyring times. The mathematical analy¬ 
sis of stochastic resonance then proceeds along the following lines. One first 
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determines the optimal tuning parameters e(T) for large T for the approxi¬ 
mating Markov chains, a rather simple task. To see that the Markov chain’s 
behavior approaches the diffusion’s in the small noise limit, spectral theory 
of the infinitesimal generator is used. Its spatial part is seen to possess a 
spectral gap between the second and third eigenvalues, and therefore the 
closeness of equilibrium distributions of the Markov chain on the one hand 
and diffusion on the other hand can be well controlled. Surprisingly, however, 
the notion of spectral power amplification is not robust for the passage from 
the Markov chain to the diffusion. Subtle dependencies on the geometrical 
fine structure of the potential function U in the minima beyond the expected 
curvature properties lead to quite unexpected counterintuitive effects. For 
example, a subtle drag away from the other well caused by the sign of the 
third derivative of U in the deep well suffices to make the spectral power 
amplification curve strictly increasing in the parameter range in which the 
approximating Markov chain has its resonance point. This dramatic devi¬ 
ation from expected behavior is due to the significance the spectral power 
amplification attributes to small intrawell fluctuations. 

Our main motivation in writing this paper was to investigate concepts 
by which on the one hand the physical intuition of optimal periodic tun¬ 
ing of random trajectories with a simple periodic input can be interpreted 
in a mathematically sound way, and which on the other hand fail to have 
this unfortunate defect of robustness. We deal with the framework of one- 
dimensional potential diffusions. The notion of quality of periodic tuning we 
shall investigate completely excludes the effect of small intrawell fluctuations 
and purely relies on the transition mechanism between domains of attrac¬ 
tion given by the potential. At the same time it generalizes the previously 
known results to potential functions which may vary periodically in time in 
a continuous, but otherwise quite general way, and whose growth at Too 



Fig. 1. Potential landscape. 
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Fig. 2. Definition of a^ 1 . 


may just be linear. More precisely, we study diffusion processes driven by a 
Brownian motion of intensity e given by the stochastic differential equation 
(SDE) 

dX t = --^u(J^,X t ^dt + VtedW t , t> 0. 

The underlying potential landscape (see Figure 1) is described by a function 
U(t, x), t > 0, x G M, which is periodic in time with period 1, and its temporal 
variation, by the rescaling with very large T, acts on the diffusion at a very 
small frequency. U is supposed to have exactly two wells located at ±1, 
separated by a saddle at 0. The depth (measured in positive quantities) 
of U(t,-) at ±1 is given by the 1-periodic depth functions D±\(t) which 
are assumed to never fall below zero. We shall throughout look at time 
scales T = exp(^), for which the Kramers-Eyring formula indicates that 
transitions, for example, from the domain of attraction of —1 to the domain 
of attraction of 1 will occur as soon as D_i becomes less than pt, that is, at 
time (see Figure 2) 

a^ 1 = inf{t > 0 : D±i(t) < p,}. 

This triggers periodic behavior of the diffusion trajectories on long time 
scales. The modern theory of metastability in time-homogeneous diffusion 
processes complements the fundamental large deviations’ theory presented 
by [8] to produce the exponential decay rates of transition probabilities be¬ 
tween different domains of attraction of a potential landscape together with 
very sharp multiplicative error estimates, uniformly on compacts in system 
parameters. Their sharpest forms are presented in some papers by Bovier, 
Eckhoff, Gayrard and Klein [3, 4], improving Day’s previous results obtained 
in [5, 6]. They are derived from deep relationships of large deviations’ the¬ 
ory with the spectral and capacity theory of the infinitesimal generator. We 
shall make use of this powerful machinery to obtain very precise estimates 
of the exponential tails of the laws of the transition times between domains 
of attraction. In fact, we have to extend the estimates by Bovier, Gayrard 
and Klein [4] to the framework of time-inhomogeneous diffusions. In the 
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underlying one-dimensional situation, this is roughly achieved by freezing 
the time dependence of the potential on small time intervals to define lower 
and upper bound time-homogeneous potentials not differing very much from 
the original one. Consequently, comparison theorems are used to control the 
transition behavior from above and below by the behavior of the correspond¬ 
ing time-homogeneous diffusions. This allows very precise estimates on the 
probabilities with which the diffusion at time scale T = exp(-) transits from 
the domain of attraction of — 1 to the domain of attraction of 1 or vice versa 
within time windows [(a^ 1 — h)T , (a^ 1 + h)T] for small h > 0. If t x {X) de¬ 
notes the transit time to i, it is shown to be given by 

lirri e ln(l — M(e,/i)) = max{/i — Z?j(a* — h)}, 
e —>0 i=± 1 M 

with 

M{e,n) = min P»(r_i(X) G [{a 1 - h)T, (a 1 + h)T}), e> 0,/r G I R , 

l=±l 

and where Ir is the resonance interval (Figure 3), that is, the set of scale 
parameters for which trivial or chaotic transition behavior of the trajectories 
is excluded (Figure 4). 

The stated convergence is uniform in fi on compact subsets of Ir. This 
allows us to take as our measure of periodic tuning, compute the 

scale fao(h) for which the transition rate is optimal, and define the stochas¬ 
tic resonance point as the eventually existing limit of po(h) as h—* 0. This 
notion of stochastic resonance is strongly related to the notions of periodic 




Fig. 4. Chaotic and trivial transition behavior of the trajectories. 
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Fig. 5. Length distribution of the interspike intervals in a simplified model (two-state 
Markov chain). 

tuning based on interspike intervals (see Figure 5 and [11]), which describe 
the probability distribution for transitions as functions of time with expo¬ 
nentially decaying spikes near the integer multiples of the forcing periods. 
As opposed to the physics notions based on spectral decomposition of the 
statistics of the solution trajectories investigated in [13] or [18] it has the 
big advantage of being robust for the passage from the diffusion to the two- 
state Markov chain reducing its behavior to the features of pure transitions 
between the two domains of attractions of metasable points. 

Here is an outline of the organization of the material in the paper. Sec¬ 
tion 1 presents a review of results from the asymptotic theory of time- 
homogeneous diffusions and their metasable sets needed for our purposes 
(Theorem 1.1). In Section 2 we bring to work the tools of comparison theo¬ 
rems to deduce the sharp exponential transition rates for our time-periodic 
diffusions from the time-homogeneous results (Theorem 2.1). Section 3 is 
devoted to applying these sharp estimates to identify stochastic resonance 
points for diffusions (Theorem 3.2), compare them to their counterparts for 
the reduced Markov chains and prove robustness of our notion of resonance 
(Theorem 3.4). 

1. Exponential distribution of transition times for time-homogeneous dif¬ 
fusions. It will turn out to be crucial for our approach of periodic tuning 
to be discussed later to obtain large deviation type estimates for the expo¬ 
nential decay rate of the law of transition times uniformly in a time scale 
parameter. We shall make use of a technique of freezing time-dependent po¬ 
tentials on small subintervals of the periodicity interval [0, T] on their states 
taken at fixed times in the intervals, to be able to use known results for time- 
homogeneous diffusions. In this setting, the uniformity problem translates 
into uniformity of the convergence to exponential decay rates in compact 
subsets of the domain of attraction the diffusion starts in and in time. It is 
clear that we are led directly into large deviations’ estimates for exit time 
distributions of time-homogeneous diffusions such as presented in the pi¬ 
oneering book by Freidlin and Wentzell [8]. But for obtaining uniformity 
in space and time, one has to use sharpened versions of these estimates 
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developed later for controlling in particular the exponential errors in the 
estimates. The purpose of this section is to summarize what we shall need 
from this fine well-established theory. 

We shall refer to the most recent and advanced development of sharp es¬ 
timates for transition times presented in [3, 4], They are valid far beyond 
our modest framework, both in the multidimensional case and for any fi¬ 
nite number of local minima of the potential. Their quality comes from a 
detailed analysis of the relationship between transition times and low-lying 
eigenvalues of the spectrum of the infinitesimal generator of the diffusion. 
We shall state them in the simple one-dimensional setting given here. A 
more complex multidimensional version can also be found in [5]. For this 
purpose, suppose that Q is a purely space-dependent C 2 potential function 
(see Figure 6) possessing only —1,1 as local minima, separated by the saddle 
point 0 at which Q takes the value 0. Suppose that the curvature of Q at — 1 
is strictly positive, that is, Q"{— 1) > 0. As for ultra- or hypercontractivity 
type properties for Q , we shall assume that it has exponentially tight level 
sets; that is, there is ao > 0 such that for any a>ao there exists a constant 
C(a) such that for e < 1, 

(1) / exp(^— dz < C(a ) expf——Y 

4{y:Q(j/)>a} V £ J \ eJ 

We shall concentrate in this situation on a transition out of the domain of 
attraction of the stable point —1 for the diffusion associated with the SDE 

dY t £ = -Q'(Y t £ ) dt + V 2e dW t , 

Y o = V- 

Let C be a closed interval of the form [d, oo[ with d ^ 0. To state our aim 
in a slightly different version, we will be interested in the asymptotics of the 
entrance time of Y £ into C : 

Tq = inf(t > 0 T/gCJ. 

Then we obtain the following result (see [4] or [5]). 



Fig. 6. Potential. 
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Theorem 1.1. Let A £ denote the principal eigenvalue of the linear op¬ 
erator 

CTu = eu" — Q'u' 

with Dirichlet boundary conditions on dC. Then for every compact K C 
] — oo, 0 [ there is a constant c> 0 such that 

( 2 ) Py(T^>t) = e- XH (l+0 K (e- c ^)), 

where Ok denotes an error term which is uniform iny G K, t> 0. Moreover, 
for the asymptotic behavior of the eigenvalue A £ the following holds: 

(3) A e Ey[rp] ■—> 1 uniformly on compacts K C] — oo,0[ 

as e — > 0. 

Proof. There are two small issues which deserve comments. 

First, the uniformity over compacts in ]— oo, 0[ claimed in the main state¬ 
ments. Day [5] tackles it. But he considers only exits from bounded domains. 
Bovier, Gayrard and Klein [4] have a version for unbounded domains, but 
uniformity over compacts of the domain of attraction in which the diffu¬ 
sion starts is not explicitly proved. It is, however, hidden in their method of 
proof of Theorem 1.3 ([4], pages 30 and 31) which makes use of an eigenfunc¬ 
tion expansion. But due to regularity results on the eigenfunctions (see [4], 
pages 16-18) they must be bounded on compacts in the domain of attraction 
]— oo,0[. This implies the desired uniformity. 

The second comment concerns our assumptions on C. Translated into 
our setting, in [4] the target set C is assumed to be closed, to contain a 
neighborhood of 1 if the potential is deeper there than at — 1 , and to have a 
positive distance from the saddle 0. Since we are in a one-dimensional setting, 
we may reduce these conditions to the simple one d 7 ^ 0. If necessary, we may 
always cut out of C a small open neighborhood of 0 without changing the 
law of tc if starting from the domain of attraction of — 1 . □ 

2. Exponential transition rates between moving domains of attraction. 

We shall now consider a potential diffusion given by the one-dimensional 
SDE 

(4) dX t = ~(± t X t ')dt + VtedW t . 

The time-periodic potential U of period 1 is supposed to fulfill the fol¬ 
lowing conditions. First of all, its global rough geometry is the one of a 
double-well potential with temporally moving wells, but time-independent 
critical points. For simplicity we suppose that its local minima are given 
by ±1, and its only saddle point by 0, independently of time. So ±1 are 
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the only metastable states of the potential on the whole time axis. Out¬ 
side of 0, ±1, is supposed to be continuous in (t,x). Our main concern 
will be the asymptotics of the transition times from the domain of attrac¬ 
tion ] — oo,0[ of —1 to the domain of attraction ]0,oo[ associated with 1 
of the time-inhomogeneous diffusion in the small noise limit e —► 0. More 
precisely, we will be interested in describing the exponential transition rate 
from the domain of —1 to the domain of 1. Our potential not being time- 
homogeneous, we shall make use of comparison arguments with diffusions 
possessing time-independent potentials in order to find a careful reduction 
of the inhomogeneous exit problem to the homogeneous one, and use the 
asymptotic results stated in Theorem 1.1 in this framework. This will be 
achieved by freezing the driving force derived from the potential on small 
time intervals on the mimimal or maximal level it takes there. To be more 
precise, for each interval I C M + let 


( 5 ) 


dU 

Vj{x) = sup —— (t,x) and 
tel ox 


dU 

Ri(x) = inf — (t,x). 
tel ox 


See Figure 7. The regularity conditions valid for U imply that V and R 
are continuous functions. Moreover, V7(—1) = Ri(— 1) = 0. If I = [a, b\, we 
denote by X the solution of the SDE on 

dx\ = -Ri{x!) dt + y/2edB t , 

( 6 ) 

Xo = X aT . 

X 1 is defined in the same way, replacing Rj by Vj. These two time-homogeneous 
diffusions are used to control the time-inhomogeneous diffusion X as long 
as time runs in the interval I. In fact, we have, P- a.s., 

X\ T < x (t+a)T < x tTi te[0,b — a\. 


Of course, to make use of the asymptotic results stated in the previous 
section, we need ultra- or hypercontractivity properties for the frozen po¬ 
tentials. To formulate a hypothesis which is both not too restrictive and 
easy to handle for time-dependent potentials, let us give the following easy 
sufficient criterion for exponential tightness of levels of a time independent 
potential Q. 

Lemma 2.1. Assume that Q is a real-valued differentiable function on 
M, and that there are constants K±,K 2 > 0 such that 

Q'(x ) < —K‘) for x < —Ki, 

( 7 ) 

Q'(x) > I \2 for x > K\. 

Then Q has exponentially tight level sets. 
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Fig. 7. Definition of Vi and Ri. 


Proof. It is obviously enough to argue on M + . Due to (7), we know that 
near oo, Q is strictly increasing with inverse Q _1 . So for a > 0 big enough 
and £>0we have 

= r “P (") WF55) dt 

-^r exp {~i) dt 

£ fa 

= ^ eXP V e 

This clearly implies exponential tightness. □ 


Conditions like (7) are practical in our setting for the following reason: 
if we assume them to be satisfied uniformly for t E [0,1], it is clear that all 
frozen potential functions will inherit the property. We will therefore assume 
throughout that there exist constants K±, K 2 > 0 such that 


dU. , 

su p 

x>o dx 

<-k 2 

for x < —K\ 

mf —(t,x) 
t> 0 dx 

>k 2 

for x> K\. 


We measure periods T on the logarithmic scale ft given by T = T e = 
exp(^). The reason for this is hidden in the classical formula of Kramers- 
Eyring. It states that in the small noise limit T is the time it takes the 
diffusing particle to climb a height /x in a potential landscape U. This for¬ 
mula has the following intuitive consequences. Assume the diffusion faces an 
obstacle of constant potential height U + > ft for exiting if it diffuses on time 
scale T. Then asymptotically it never exits on this scale. On the other hand, 
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if it faces an obstacle of height E/_ < /i diffusing at time scale T, it has to 
exit immediately in the small noise limit. We introduce the depth function 
at x E M by 


D x (t) = U(t,0) — U(t,x), t> 0. 

The maximal well depths Z)_i and D\ will be of particular importance. We 
shall assume that they satisfy the assumption (see Figure 2) 

all local extrema of D\, D-\ are global and the 
(9) functions are strictly monotonous between the ex¬ 

trema. 

For \i > 0 let now 


an = inf{t > 0: D_i(t) < /i}. 

This is the same function as a~ l defined above. Here we omit the superscript 
since we always concentrate on transitions from —1 to 1. The interval 

( 10 ) I-i = 

contains all possible depths the potential minimum located at —1 takes 
during one period of time. Note that = 0 for /r strictly above the upper 
boundary of /_i, and = oo for /j, strictly below the lower boundary. If 
/i € I— i, then for times beyond a^, the barrier height between —1 and 0 
has dropped below the critical level so that on time scale T the diffusion 
exits immediately. According to this heuristic, the diffusion running on time 
scale T should exit the domain of attraction ] — oo,0[ through 0 and then 
transit to the other well immediately after TD_\(t ) drops below /i. We shall 
be interested in the exponential rate at which this happens, uniformly in 
starting points taken from a compact in ] — oo,0[. For this purpose, for a 
regular diffusion Y on R and a £ 1, we denote by r a (Y) the first hitting 
time of a by Y. 

Our first aim is to prove the precise estimate 

(11) limelnP x (n(X) < (a M - h)T) = n~ T>_i(a M - h) 

for x < 0, uniformly in /r on compact subintervals T of (10) and for 0 < h < 
inf^gr Q/f Note that given such a compact F, due to the continuity of 
we have inf^gra^ > 0. Fix T,x and h from now on. Our approach proceeds 
in essentially two steps. 


inf i (t), sup D-\(t) 
*>° t>o 
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2.1. Upper bound for the exponential exit rate. In the first step we shall 
find upper bounds for the exponential rates. For this purpose we shall 
partition the relevant time intervals [0, — h\. Fix some 5 > 0 such that 

5 < |x|. Since ^ is continuous, we may choose an equidistant partition 
0 = ro < 7 T < ■ • ■ < r n = — h of [ 0 , — h\ with mesh 7 small enough to 

ensure 


( 12 ) 


sup sup 

s,t£[rj,rj. |_i] x€[— 1,0] 


dU. , 


du , , 


<5. 


Denote Ij = rf\, Though the choice of the intervals depends 

on h and T, their number will be bounded by a universal constant. Using 
this partition, we may start our search for an upper bound by freezing the 
time-dependent potential on Ij at its value Ri j , 1 < j < n, and then com¬ 
paring the diffusion X there with X Ij . There is a little difficulty with this 
procedure. The drift coefficients Rj which govern the SDE for X Ij were 
defined by taking infima over time intervals. This operation may destroy 
their differentiability properties in the spatial variable. Therefore it may be 
necessary to compare X on the intervals Ij with smoother diffusions still 
dominating it. But this can be done at no extra cost. For each 1 < j < n we 
may choose Rj EC 1 )!!) satisfying 


Rj 5; Rij L Rj + h, 

there are rrij E] — 1 — 6, — 1 + <5[, Sj E] — 5, S[.rrij G]1 — 6 ,1 + <5[ 

/-.Q\SUch that Rj\]—oo,m, .[Uls,,m,f ^ 0; ^j\]m ,Si [U lm,' ,oo[ 0) 

[ ) Rj{rn j )> 0 , 

there are constants K±, K 2 > 0 such that Rj(x) < — K 2 for x < —K\. 
Rj(x) > K 2 for x > K\. 

Let Xi be the diffusion associated with Rj , 1 < j < n. Let us choose a par¬ 
tition x = xq < x\ < ■ ■ ■ < x n = —6 of the interval [x, —5] which will typically 
not be supposed to be equidistant. In the following inequality the diffusions 
X and X ] on Ij are compared and the Markov property is employed. We 
have 


P*(ri (X)<{a„-h)T) 

n— 1 

< ^ > x{r Xj _ 1 (X) > rj-iT,T Xj (X) < rjT) 

3 = 1 

+ Pa;(Tx„_ 1 (-Y) > r n -\T,Ti(X) < r n T) 

n— 1 

< E Vx[{T Xj (X. +rj _ lT ) < 7 T} n {T Xj _i{X) > rj -\T}\ 

3 = 1 
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(14) 

+ P x [{n(X +r . ri _ lT ) < 7 T} n (X) > r n -iT}] 

n— 1 

< ^E x [P A - r ._ iT (r Xn (X J ) 
i=i 

+ E x [Px rn _ lT (ri(X n ) < 7 r ) 1 {T BB _ 1 (A-)>r„_ 1 r}] 

n—1 

< E ? -7 :(-d(- V '') < 7^) +Px„_ 1 (r 1 (X n ) < 7 T). 

1=1 

Let us now fix 1 < j < n and continue estimating the terms Px -_i(r x - (X- 7 ) < 
7 T) and P^.! (ri(X n ) < 7 T) individually. For this purpose we apply Theo¬ 
rem 1.1 for Q = Rj,d = Xj to obtain that 

P Xj . 1 (T Xj (X^)<'yT) 

< 1 - e" A ) 7T (l - e ~ c/£ ) < 1 - e" A ) 7T + e~ c / £ , 

(15) 

Px„_i ( T i(X n ) < 7 T) 

< 1 - e" A " 7T (l - e~ c / £ ) < 1 - e" A » 7T + e" c/£ , 

uniformly in T, hence uniformly in /r. Here A^ denotes the principal eigen¬ 
value of the operator C J defined by 

C £ jU = eu" — 

with Dirichlet boundary conditions at Xj. We now come to the crucial part 
of the derivation of an upper estimate. We shall use precise asymptotics of 
the eigenvalues A^. 


Lemma 2.2. There exists C > 0 such that for 1 < j < n, 


(16) 


lim e In A^ — [ U(rj,Xj ) — U(rj , 


D] 


<CX 


Proof. Fix l<j<n. Define the pseudopotential corresponding to the 
drift Rj by 

Vj{x,z) =inf|i^ ($, + -y(0s)) ds,4> 0 = x,(j)t = z,t>0 

where (j) stands for absolutely continuous functions defined on the time in¬ 
terval [0,t\. Since due to our assumptions m J is the only local minimum of 
the potential corresponding to Rj on ] — oo,0[, the sharpened form of the 
exit time theorem of Freidlin-Wentzell (see [4], Theorem 1.1) implies that 
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Let us estimate the pseudopotential. We have 
Vj(rrij,Xj)= I Rj(9)d0 


[ x i dU 
/-1 Ox 


fa,9)d0 + 


BU 

R i(°)-^(r j ,0))d0 


LJI^ e)de - 


Continuity of in (t, x ) entails the existence of Ci < 0 such that 


To estimate the second remainder term, recall that the mesh 7 was chosen 
to produce at most 5 as modulus of continuity of ^ [see (12)], and that Rj 
is also at most a distance 6 away [see (13)]. We therefore obtain 

< 26. 


Hence 

(17) 


R jm 


lim e In A j - [Ufa,Xj) - Ufa,- 1)] 


3 

The asserted asymptotic result follows. □ 


< (2 + C\)5. 


As a consequence we obtain an upper bound for the exponential conver¬ 
gence rate for the exit time from the domain of attraction of the potential 
well at — 1 . 


Proposition 2.1. Let x<0 and let T be a compact subset of (10). 
Then there exists 0 < ho < inf Me r such that for h< fa, 

(18) lim sup e In (ri (X) < (a M — h)T) < fi — H_i(a At — h) 

£—>0 

uniformly for fx E T . 

Proof. According to what has been proved, there are constants Eq > 0 
and K > 0 such that for e < eq, /x E T, 

P*(ri(X) < (a M - h)T) 

n— 1 

< E v Xj faT Xj (x 3 ) < tT) + P a ,„_ 1 (n(x") < 7 t) 

3 =1 

n 

< K^T E + ne~ c / £ . 
i=i 


(19) 
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Taking logarithms on both sides, multiplying by e and using the equation 
lirrieln(/(e) + g{s)) = max limeln(/(e)), limeln(g((e)) 

[ e —>0 

for two positive functions / and g, we may apply Lemma 2.2 to get 

limsupelnP x (ri(X) < (a M — h)T) 

£—>0 

< [max[/r — [U ( rj,Xj ) — U ( rj , —1)]: 1 < j < n] + CS\ V (—c) 

< [max[g — [U (rj,xi) — U{rj, — 1 )]: 1 < j < n] + C5\ V (— c). 

Recalling the definition of a^ and that x\ < 0 is arbitrary, we may conclude 

limsupelnP x (n(X) < (a M — h)T ) 

e —>0 

< [max[/j — D_i(rj): 1 < j < n] + C5] V (— c) 

< [(n - -D-i (a M - h)) + C(5] V (-c) 

uniformly for /igT. Now choose ho > 0 small enough so that for h < ho we 
have 

inf (fi — — h)) > —c. 

/ier 

Finally, since 5 is arbitrary, we may let 5 tend to zero. This way we obtain 
the desired upper bound for the exponential rate. □ 


2.2. Lower bound for the exponential exit rate. In the second step of our 
approach, we shall establish lower bounds for the exponential rates at which 
the diffusion exits from the basin of attraction of —1. Let us first prove an 
auxiliary result. It states that the probability of exiting the interval [1,0] 
via l is exponentially small with exponential order increasing in \l\, due to 
hypercontractivity. Recall the constants I\± and K 2 from (8). 


Lemma 2.3. There exist positive constants C and £q such that for e < £0 
and l < x A —K\ and g > 0 we have 


( 20 ) 


IPx(t;(X) < T) < — exp( 


C {2K 2 (l-xA(-K 1 )) + n 


Proof. We give arguments for the case x < —K\, the other case being 
easier. Recalling that by (8) the gradient of U is bounded below by —K 2 on 
R+ x] — 00, — Ki] j we may compare the diffusion X with the diffusion Z on 
the interval ] — 00, x] reflected at x with constant drift equal to K 2 . It can 
be given by the SDE 

dZt = K 2 dt + dL t + y/2e dWt , 
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where Zq = x and L is an increasing process satisfying Jq(Z s — x) dL s = 0, 
t > 0. See, for example, [19]. The comparison clearly yields 

(21) P*(7l(X)<T)<P x (^(Z)<T). 


By Chebyshev’s inequality, 

(22) ¥ x (ri(Z) <T) < eE x 


exp( --n{Z) 


Let ip{y) = E ?/ [exp(—^r;(Z))], l <y <x. Our task consists in an estimation 
of <p(x). According to the Feynman-Kac and Dynkin formulae ip solves the 
boundary value problem 


Eip" + K 2 ip' - = 0 

ip' {x) = 0, ip(l) = L 


on ]l, x [, 


The eigenvalues of the differential equation are determined by the equation 

eA 2 + K 2 X — — = 0, 

hence by A^ = ^[— K 2 ± yjK‘j + ^]. Taking the boundary conditions into 
account leads to the equation 

\+ \+x+\-y _ \~x+X+y 

= ^ +e \+ x+ \-l _ y- e A-x+A +1 ’ y£[li x \- 

Neglecting the second term in the denominator of the fraction, we obtain 


ip{x) < 


\J KI + 4e/T 

E\ + e x ~( l ~ x ') 


Now for e small enough, A + > and A < —^. Therefore, for e small 
enough there exists a constant Cq > 0 independent of e and T such that 


<p(x ) < CoK^ — exp 


\K 2 {l-x)^ 

CoKl 

=i 

+ 

TT 

1 

'■■-o 

£ 

— -exp 

£ 

£ 


This implies the desired inequality (20). □ 


We shall continue to use the partition (Ij : 1 < j < n) of the interval 
[0, an — h] of the preceding section. This time, we shall compare with ho¬ 
mogeneous diffusions by freezing the potential derivative at an upper level, 
which results in working with the drifts V Ij and the diffusions 2L Ij ■ Again, 
these drifts may fail to possess the regularity properties required to apply 
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Theorem 1.1. For this reason we may choose smoothed versions Vj S C 1 (R) 
of the potentials with corresponding diffusion processes Y J satisfying 

Vj > V Ij >Vj + 5, 

there are rrij € ] — 1 — <5, — 1 + <5[, Sj € ] — 8, <5[, nrij £ ]1 — 5, 1 + <5[ 

that Vj |]_oo,m , [U]sj,mj[ ^ 0? Vj llm ■■s-.lUlTO-.-.oof 

{ 0 , 

there are constants K i. K 2 > 0 such that Vj(x) < —K 2 for x < —Ad, 
Vj(x) > K 2 for x> K\. 

To deduce a lower estimate, we shall compare X via X In with Y n on the 
interval I n in the scale T. We may write for l < x 

^ x (T 1 (X)<(a li -h)T) 

> P x (T!(X) < r n T , n (X) A n(X) > r n _\T) 

(24) 

>E,(l {n( x)Au W >, n _iT}^ rn _ lT (ri(X +rri _ lT ) <7 T)) 

> P i(n(Y n ) < 7 T) x P x (n(X) A n (X) > r n -\T). 

As a consequence of /r — Z)_i((a At — /i)T) < 0 and the arguments presented 
in Section 2.1, we note that uniformly on our compact set T 

limP x (n(X) < r n -iT) = 0. 

£—>0 

This clearly implies that there is £0 > 0 an d a constant C > ^ such that for 

£<£o, 

(25) ^ x {Ti{X)>r n -iT)>C. 

Moreover, by Lemma 2.3, for l small enough, there exists £\ > 0 such that 
for e < El, 

P*(ri(X) > r n -iT) > C. 

Hence for e small enough we have 

P x (n(A) A n(X) > r n -\T) > c - i > 0. 

It therefore remains to find lower bounds for P/(ri(T n ) < 7T). 

We may now apply the same arguments as those developed in Section 2.1. 
We just have to use Lemma 2.2 for the eigenvalues of the operator 

C e u = eu" — V n u! 

with Dirichlet boundary conditions at 0 in the sense of lower bounds uni¬ 
formly on the compact T. As a consequence we obviously obtain, with a 
constant C > 0 independent of T, 

liminf elnP;(ri(T n ) < 7T) > (/z — H_i(a /i — h) — CS) V (—c), 
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uniformly for /iGT. Let us now choose ho > 0 and do > 0 small enough such 
that for h < ho,5 < do we have fx — D_i((a M — h)T) — C5 > —c. Since 5 is 
arbitrary, we obtain 

liminfelnP z (n(y n ) < yT) > (/r - L»_i((a M - h)T)), 

uniformly for /i£T. Recalling (24) and (25), we finally obtain 

liminfelnP x (n(X) < (a M - h)T) >n~ L>_i((a M - h)T), 

uniformly on T. With this result we have established the desired lower bound 
for the exponential exit rate. 

Proposition 2.2. Let x < 0, and let T be a compact subset of (10). 
Then there exists 0 < ho < inf p _ e r such that for h < ho, 

(26) limelnP x (n(X) < (a M - h)T) > //- D_i(a M - h) 

uniformly for /i £ T. 

2.3. The exponential smallness of the rate of too long transitions. Having 
proved (11) in the preceding two propositions, the second aim of this section 
is to show that the exponential rate at which T\ (X) exceeds (a M + h)T 
is arbitrarily small. In fact, we shall make precise that the rate at which 
transitions happen which take at least as long as (a^ + h)T vanishes to all 
exponential orders, for h > 0 arbitrary. 


Proposition 2.3. Let x < 0, and let T be a compact subset of (10) not 
containing D_ i(0). Then there exists ho > 0 such that for all 0 < h < ho and 
p, E T roe have 

limsups:lnP x (Ti(X) > (a^ + h)T) = — oo. 

£—>0 


Proof. Let 5 > 0 and h > 0 be given. Let T be a compact subset of 
(10). First, for l < x we may write 

(27) P*(ri(X) > ( fl/1 + h)T) < PsMX) A n (X) > (a M + h)T) 

(28) + Pi(r(X) < (a^ + h)T). 

To estimate the second term on the right-hand side of (27), we employ 
Lemma 2.3. In fact, for l < x A (— K{) we have 

lirri sup £ In P x (r; (X) < (a M + h)T) < 2 X 2(1 — x A (— K\)) + sup^r. 
e^O r 

Therefore 


lim 

l—> — 00 


limsupelnP x (r;(X) < (a M + h)T) 
. £—>0 


= — OO. 
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It therefore remains to estimate the first term on the right-hand side of 
(27) for l small but fixed. Let 0 = ro < r\ < ■ ■ ■ < r n = + h be an equidis¬ 
tant partition of the interval [0, + h] of mesh 7 < 4 and denote Ij = 

, i'j], 1 < j < n. Then we have 

W x (T l (X)/\Ti{X) >{<1^ + 11)^ 

= F x (t 0 (X) Ati(X) > r n -iT,Ti(X. +rn _ lT ) A Ti(X. +rn _ lT ) > 7 T) 

< K x(^{T 0 (X)^Ti(X)>r n ^ 1 T}^ , X rri _ 1 T( T l{2L In ) > 7 T)) 

(29) 

< max r P 2 / (ri(^)> 7 T) 

= Pi(n(X / ")> 7 T) 

<P J (r 1 (X n )> 7 T). 


Here, we compare the inhomogeneous diffusion X on I n with the time- 
homogeneous one X In corresponding to drift Rj n and finally with X n subject 
to drift R n to be described below, and we use monotonicity of 

y^Py(ri(X J ")>7T). 


We assume 7 to be small enough to ensure 


(30) 


sup sup 
s,t£[r n -i,r n ] xe[Z,0] 


dU_ 

dx 


(t,x) 


dU_ 

dx 


(s,x) 


<5. 


We may choose the drift R n to satisfy 


Rn Rl n < Rn + 8, 

there are m n £ ] — 1 — 8, — 1 + d[, s n £ ] — 8,8[,m n £ ]1 — 5 ,1 + <5[ 

/oI \Such, that Rn\]— [uls„.m„,[ ^ 0, l]m n ,Sn[U]m„,oo[ 0) 

1 X(™j>o, 

there are constants iTi, A '2 > 0 such that R n (x) < —K 2 for x < —K\ ■ 

Rn (x) > K 2 for x > K\. 

We are ready to apply Theorem 1.1, this time for Q = R n , d= —5 to obtain 
that 

(32) Pz(ri(X n ) > 7 T) < e~ x ^ T (l + e ~ c/e ) 


uniformly on T. Here \ £ n stands for the principal eigenvalue of the operator 
C £ n defined by 

C £ n u = eu" — R n v! 


with Dirichlet boundary conditions at —(5. The asymptotic properties of A^ 
can be deduced in a similar way to Lemma 2.2. We estimate the pseudopo¬ 
tential corresponding to R n taking (30) into account. We obtain that there 
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exists C > 0 such that 


(33) 


limeA^T- (/i-D_i(r n )) 

£—>0 


< C5 


uniformly on F. Now recall that due to the choice of 7 , we have < r n , 

hence 


/r- -D-i(r n ) > 0, /iGT, 
and by compactness of T even, 

inf [fi - D-i(r n )} > 0. 

This in turn implies that 

lim X £ T = 00 

£—>0 n 

uniformly on T. But due to (32) we are allowed to conclude 

lim P/(Ti(X n ) > 7 T) = —00 
e —>0 

uniformly in T. According to (29), this completes the proof. □ 

Finally, we may summarize the results of Sections 2.1 and 2 . 2 , and state 
the main result on asymptotic exponential decay rates of transition proba¬ 
bilities. 

Theorem 2.1. Let x < 0, and let T be a compact subset of (10). For 
e > 0, p>0 let T = exp(^). Then there exists h$ > 0 such that for h < ho, 

(34) linrelnP x (ri(X) £ [(a M - h)T, (a M + h)T]) = /r -D_i(a M - h) 
uniformly for /i£T. 

Proof. This is an immediate consequence of Propositions 2.1-2.3. □ 

3. Stochastic resonance in a double-well potential. Let us now turn to 
the main subject of this paper, a characterization of the notion of stochastic 
resonance. Let us recall that we look for a characterization of the concept 
of optimal periodic tuning which is extensively studied in the physics lit¬ 
erature by notions such as the signal-to-noise ratio or the spectral power 
amplification (see [18]). Let us also remark that this concept implicitly uses 
and refines the concept of stochastic resonance studied by Freidlin [7] which 
paraphrases the ability of periodically perturbed stochastic systems to follow 
the periodic excitation in the small noise limit, and exhibit quasi-periodic 
motion. In more mathematical terms and the notation introduced before, we 
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aim at choosing the noise intensity parameter e such that in the large period 
limit T —► oo the diffusion trajectories follow the periodic excitation of the 
system hidden in U in an optimal way to be made precise. In Section 3.1 we 
shall show that a quality measure of goodness of periodic tuning is given by 
the exponential rate at which the first transition to the other well happens 
within a fixed interval around a^T. In Section 3.2 we establish robustness 
of this notion of quality: we show that in the small noise limit the diffusion 
and its reduced model, a Markov chain living on a two-point state space, 
have the same resonance pattern. 


3.1. Transition probabilities as a measure of quality. The local extrema 
of the depth functions D±\ of U are supposed to be global, and D±\ is 
strictly increasing between its extrema. Recall that we work with exponential 
time scales /u related to the natural time T by the equation T = exp(^). In 
this section, we have to work with scale functions depending on the starting 
well and eventually on arbitrary starting times. So we let 

a /j,( s ) = hif{i > s : Di(t ) < //}, i = ±1, n > 0. 

The relevant time scales /i will be chosen from the intervals 


h = 


inf AO), sup Di(t) , 
t>o t >o 


i = ± 1. 


Our aim is to observe periodic behavior of the diffusion. This will in principle 
mean that the process can travel from one well to the other and back on the 
time scales in which we let the diffusion run, but not instantaneously. So, 
on the one hand, we have to work on time scales on which it does not get 
stuck in one of the wells of the potential. On the other hand, the time scales 
we are concentrating on should also not allow for chaotic behavior , that is, 
immediate re-bouncing after changing the well. 

To make these conditions mathematically precise, recall that transitions 
become possible as soon as the potential barrier D ±i becomes smaller than 
the time scale parameter //. Hence if /i > infj>o Di(t), there is a time range 
during which the diffusion can leave the well centered at i. To not get stuck 
in one of them, the diffusion has to be able to leave both. This is guaranteed 
if 


(35) /i > max inf D* (t). 

i=± 11>0 

To avoid immediate re-bouncing, we have to assure that the diffusion 
cannot leave the domain of attraction of —i at the moment it reaches it, 
coming from i. Suppose we consider the dynamics after time s > 0, and the 
diffusion is near i at that time. Its first transition to the well at —i occurs 
at time a^(s)T, and it stays there for at least a little while if -D_j(a^(s)) 
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is bigger than /r. This is equivalent to stating that for all s > 0 there exists 
5 > 0 such that on [a^(s),a^(s) + 5] we have // < D-i. But for t shortly 
after a* (s), we always have Di(t) < ft by the very definition of a 1 ^. Hence 
our condition becomes equivalent to the following: for all s > 0 there exists 
5 > 0 such that on [a^(s), a^(s) + 5] we have ft < rnax 7 ; = -|-i D{. This in turn 
is more elegantly expressed by 


(36) ft < inf max Di ( t ). 

t>0 i=±l 

See Figure 8. 

We may summarize our search for an appropriate set of scale parameters ft 
for which periodicity in the diffusion behavior will occur. We call this set the 
“resonance interval” to indicate that we have to look for the scale of optimal 
periodicity, the resonance scale, in this interval. See [11] for the definition 
of the corresponding interval in the case of two-state Markov chains. The 
interval 


Ir 


max inf DAt ), inf max Dj (t) 
i=±it>o v J ’t>0i=±l n ’ 


is called the resonance interval (see Figure 3). Let us pause for a moment 
at this point to compare our approach with Freidlin’s [7] understanding 
of stochastic resonance by quasi-deterministic motion. In Freidlin’s terms, 
stochastic resonance is given if the parameter /r exceeds the lower boundary 
of our resonance interval. Our concept of resonance stipulates to look for an 
optimal /i in the resonance interval at which in a sense to be made precise 
the quality of periodic tuning is optimal. 


Let us now come to the discussion of the quality of periodic response of 
the stochastic system given by the diffusion, in dependence on the noise 
parameter e and the time scale parameter r which according to the remarks 
made above has to be chosen in the resonance interval. To simplify things a 
little, let us assume that the depth functions are related by a phase 4> ]0,1[, 

that is, 

D-i(t) = Di(t + (f>), t> 0. 


D-x(t) D,(i) 



Fig. 8. Depth functions in phase. 
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Moreover, we assume that the diffusion starts in —1. There are many ways 
to describe optimality of periodic tuning. Imkeller and Pavlyukevich [12, 13] 
consider different measures of quality such as the spectral power amplifi¬ 
cation, the energy, the energy-to-noise ratio, the out-of-phase measure, the 
entropy and the relative entropy. The detailed study of the physicists’ fa¬ 
vorite measure, the spectral power amplification, based on the energy of the 
spectral component of the mean trajectory in equilibrium corresponding to 
the forcing frequency ffi, shows one surprising defect: it is not robust as one 
passes from the diffusion to a reduced model described by a two-state Markov 
chain jumping with rates corresponding to the transition rates between the 
metasable states ±1 of the diffusion given by the potential minima. In fact, 
while the Markov chain’s spectral power coefficient shows a pronounced peak 
for p near an average well depth, the overwhelming influence of the diffu¬ 
sion’s fluctuations in small neighborhoods of the potential wells, discovering 
very subtle details of the potential’s geometry there, destroys this picture 
completely. Here we propose a notion of quality of periodic tuning which is 
based on the pure transition mechanism of the system between the domains 
of attraction of the double-well potential. Generalizing the approach of a 
study of optimal tuning for two-state Markov chain models (see [11]), we 
measure the quality of tuning by computing for varying time scale param¬ 
eters p the probability that, starting in i, the diffusion is transferred to — i 
within the time window [(a^ — h)T, (a^ + h)T} of width 2 hT. To find the 
stochastic resonance point for large T we have to maximize this measure of 
quality in p £ I R . The probability for transition within this window will be 
computed by the estimates of the preceding section. Uniformity of conver¬ 
gence to the exponential rates will enable us to maximize in p. Let us now 
make these ideas precise. 

To make sure that the transition window makes sense at least for small h, 
we have to suppose that > 0,i = ±1 for p € Jr. This will be guaranteed 
if 

(37) -Dj(O) > inf maxllj(f), i = ± 1. 

t>0 i=±l 

If this is not granted from the beginning, it suffices to start the diffusion a 
little later, in order to be sure that (37) is satisfied. Under (37), we call 

(38) M(e,p) = mmPfiT-fiX) £[(a l -h)T,(a l +h)T]), £>0,p<El R , 

l=±l ^ ^ 

the transition probability for a time window of width h. 

We are prepared to state our main resonance result. 

Theorem 3.1. Let T be a compact subset of I R , and let ho > 0 be given 
according to Theorem 2.1. Then 

(39) limeln(l — M{e,p)) = ma x{p — Dfiaf — h)} 

e^O i=± 1 M 
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uniformly for /i£T. 

Proof. This proposition is an obvious consequence of Propositions 2.1- 


It is clear that for h small the eventually existing global minimizer plr(K) 
of 

Ir3 /in max{/j - Dfa 1 - h)} 

l=±l ^ 

is a good candidate for our resonance point. But it still depends on h. To 
get rid of this dependence, we shall consider the limit of nn(h) as h —> 0. 

Definition 3.1. Suppose that 

Ik9/iw max{/z — Di(a l u — h)} 
i=±l ^ 

possesses a global minimum pt r(S ). Suppose further that 

Hr = lirn p R (h) 
h —>0 

exists in Ir. We call the stochastic resonance point of the diffusion X 
with time-periodic potential U. 

We shall now show that the stochastic resonance point exists if one of the 
depth functions, and thus both, due to the phase lag, has a unique point of 
maximal decrease on the interval where it is strictly decreasing. See Figure 9. 


Theorem 3.2. Suppose that D\ is twice continuously differentiable and 
has its global maximum at t\, and its global minimum at t 2 , where t\ <t 2 - 
Suppose further that there is a unique point t\ < s < t 2 such that Di|i tl S r 
is strictly concave, and Di\\ s ,t 2 \ strictly convex. Then p.R = D\{s) is the 
stochastic resonance point. 



Fig. 9. Point of maximal decrease. 
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Proof. First of all, note that there is g] 0,1[ such that D\ = D- 1 (- + 
ip). As a consequence of this, 

max{/r - Di(a l - h)} = {/i - D^a 1 - h)}. 

Z=±l 

Write = a* and recall that on the interval of decrease of D\, a M = Dj” 1 (/r). 
Therefore, the differentiability assumption yields 

1 = D'i{aft - h ) • a'n = L»i(a M - h) ■ j . 

Our hypothesis concerning convexity and concavity of D\ essentially means 
that D'{(s) = 0, and D" |]i ljS [ < 0, D”\j s t2 ^ > 0; in other words, that /i i—> 
D'^a^) has a local maximum at a^ = s. Hence for h small there exists a 
unique point a^{h) such that 

D[(a^h) - h) = D[(a^h)) 

and 


lima /t (7i) = s. 

h—>0 

To show that a^{h) corresponds to a minimum of the function 

D 1 (a fl - h)\, 

we take the second derivative of this function at a^(/i), which is given by 

D[(a^h) - h)D'{(a»(h)) - D'{(a^h) - h)D[(a^(h)) 

D[Mh)) 

But D[(a^(h)), D'^a^Qi) — h) < 0, whereas D'[(a^(h) — h) > 0, D'[(a^h)) < 
0. This clearly implies that a M (h) corresponds to a minimum of the func¬ 
tion. But by definition, as h —► 0, a^{h) —► s. Therefore, finally, D\(s) is the 
stochastic resonance point. □ 

To illustrate our results, we next discuss an example. 

Example. Let us consider the double-well potential 

U (t, x ) = Y~ cos { 27r (* ~\+ V'Sgn^)^) j (y - y 

with T = exp(^) and ip G [0, |[. See Figures 10 and 11. 

U satisfies all the assumptions required for potentials above, in particular 
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Fig. 11. Level sets of the potential. 
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—1 and 1 are the metasable states of the potential, and 0 is the saddle point. 
We can then compute the barrier height of the two wells. For i S {—1,1}, 

,,24/ tt\ 

Di(t) = -- i— cos( 2-7rt + *27n/> - - j. 

Let us note that Di(t) = D_i(t + 2 tp + ^). Since ^ S [0, 4[, we are in the 
phase case with <f> = + 1/2. The resonance interval is then given by 


Ir 


224 (n r 

—,-cos — 

5 3 15 V 2 



In the symmetric case, that is, if -0 = 0, we obtain Ir =]2/5,2/3[. Let us 
now compute the optimal tuning scale applying Theorem 2.1. We obtain 

-l 1 /15//i 1\\ , 1 

=^“ rc ” s (,Y(,2-3jj +, '’ + 4' 

Hence, for h > 0 small enough, 

=V-D-i(a ~ 1 - h ) 

="-1 - te ™( arccos {f (j - 5 )^ ^ 

/ 2\ 4 

= (//——j(l — cos 27r/i) — — sin(27rh) 



Let us recall that F does not depend on the phase which implies 

r - D.^a- 1 - h)T ) = n- - h)T). 

Hence, to obtain optimal tuning, it suffices to compute the minimum of F 
for fi E Ir. Differentiating F, we obtain 


F'{fj) = 1 — cos( 27r/i) H—- sin(27rh) 

W y 1 2 y ’ -^/l — (15/x/4 — 15/6) a 

Hence F attains its minimum for 


R R {h) = | - -cos(2vrh) 

and 


rr= lim n R (h) = §. 

h—> 0 

Thus we obtain that rr is the stochastic resonance point if rr G Ir, that 
is, if the phase is near to ^, that is, if is close to 0. In the other case, the 
optimal tuning rate on every interval [a, b] C Ir is given by the upper bound 
b. 
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3.2. The robustness of stochastic resonance based on transition windows. 
In the small noise limit e —> 0, it seems reasonable to assume that the peri¬ 
odicity properties of the diffusion trajectories caused by the periodic forcing 
due to the potential, are essentially captured by a simpler, reduced stochas¬ 
tic process: a continuous-time Markov chain which just jumps between two 
states representing the bottoms of the wells of the double-well potential at 
rates corresponding to the transition mechanism of the diffusion. This is 
just the reduction idea ubiquitous in the physics literature, and explained, 
for example, in [15]. In [13] it is found that this idea may conflict with the 
intrawell fluctuations of the diffusion if the quality of periodic tuning is mea¬ 
sured by concepts using spectral decompositions of the trajectories. We shall 
now show that in the small noise limit both models, diffusion and Markov 
chain, produce the same resonance picture, if quality of periodic tuning is 
measured by transition rates as discussed in Section 3.1. 

We first have to describe the reduced model. Let U be a time-dependent 
potential function generating the potential diffusions of the preceding sec¬ 
tion. Recall that the depth functions of the potential minima satisfy D\(t) = 
D-i(t + 4>), t > 0, with phase shift 4> E]0,1[. So, let us consider a time- 
continuous Markov chain {Y t ,t > 0} taking values in the state space { — 1,1} 
with initial data Yq = — 1. Suppose the infinitesimal generator is given by 



where ip(t) = (p(t + (j>), t > 0, and <p is a 1-periodic function describing a rate 
which just produces the transition dynamics of the diffusion between the 
potential minima ±1, that is, 

(40) y?(f)=exp^- t> 0. 

Note that by choice of <f>, 

(41) i/j(t) = exp^— t > 0. 


Transition probabilities for the Markov chain thus defined are easily com¬ 
puted. See ([11], Section 2). For example, the probability density of the first 
transition time Uj(T) is given by 


p(t) = <p(t) exp^— J ip(s)ds 
q(t) = ip(t + 4 >) exp J (p(s + (j>) ds 


if i = — 1, 


if i = 1, 
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t> 0. Equation (42) can be used to obtain results on exponential rates of 
the transition times if starting from —i, i = ±1. We summarize them 

and apply them to the measure of quality of periodic tuning in case (37): 

(42) N(s,n) = minPj(cr_j(y) € [(a? -h)T, {a 1 +h)T]), 

l=±l ^ ^ 

which is called transition probability for a time window of width h for the 
Markov chain. 

Here is the asymptotic result obtained from a slight modification of Theo¬ 
rems 3 and 4 of [11] which consists of allowing more general depth functions 
than the sinusoidal ones used there and requires just the same proof. 

Theorem 3.3. LetT be a compact subset of Ir and let ho <sup(a / J 1 ,T/2 — 
aff 1 ). Then for 0 < h < ho, 

(43) lirri £ ln(l — N(e,p,)) = max{/i — D t (a) l — h)} 

e —>0 i=± 1 ^ 

uniformly for /i£T. 

It is clear from Theorem 3.3 that the reduced Markov chain Y and the 
diffusion process X have exactly the same resonance behavior. Of course, 
we may define the stochastic resonance point for Y just as we did for X. So 
the following final robustness result holds true. 

Theorem 3.4. The resonance points of X with periodic potential U and 
of Y with exponential transition rate functions D ±i coincide. 
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